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Abstract The X-ray spectral and timing properties of 
ultraluminous X-ray sources (ULXs) have many simi- 
larities with the very high state of stellar-mass black 
holes (power-law dominated, at accretion rates greater 
than the Eddington rate). On the other hand, their 
cool disk components, large characteristic inner-disk 
radii and low characteristic timescales have been in- 
terpreted as evidence of black hole masses ~ IOOOMq 
(intermediate-mass black holes). Here we re-examine 
the physical interpretation of the cool disk model, in 
the context of accretion states of stellar-mass black 
holes. In particular, XTE J1550— 564 can be consid- 
ered the missing link between ULXs and stellar-mass 
black holes, because it exhibits a high-accretion-rate, 
low-disk-temperature state (ultraluminous branch). 
On the ultraluminous branch, the accretion rate is 
positively correlated with the disk truncation radius 
and the bolometric disk luminosity, while it is anti- 
correlated with the peak temperature and the frequency 
of quasi-periodic-oscillations. Two prototypical ULXs 
(NGC 1313 X-1 and X-2) also seem to move along that 
branch. We use a phenomenological model to show how 
the different range of spectral and timing parameters 
found in the two classes of accreting black holes de- 
pends on both their masses and accretion rates. We sug- 
gest that ULXs are consistent with black hole masses 
~ 5O-IOOM0, moderately inefficiently accreting at « 20 
times Eddington. 
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1 Introduction 

Black holes (BHs) are very simple systems, entirely 
characterized by mass and spin. However, accretion 
onto BHs is a much more complex process. Observa- 
tionally, accreting stellar-mass BHs are found in vari- 
ous transient states, with different spectral and timing 
properties, and a different balance between mass-energy 
advection, radiation and outflows. For BHs in a binary 
system, accreting from a Roche-lobe-filling donor star, 
the main parameter that determines the accretion state 
is thought to be the mass accretion rate, or more pre- 
cisely the mass inflow rate M at the outer edge of the 
accretion disk. This can be larger than the accretion 
rate through the BH horizon, if some of the gas is lost 
in outflows. The inflow rate depends on the nature and 
orbital separation of the companion star. It is conve- 
nient to rescale this rate as a dimensionless quantity, 
m = M/M^dd, where Mfidd = iEdd/(0.1c^) is the ac- 
cretion rate that would be required to produce a lumi- 
nosity L — L^dd — 1-3 X 10'^* erg s~^ with a standard 
radiative efficiency = 0.1. 

BHs accreting at a rate m < a few 10~^ are generally 
found in the low/hard state (Remillard and McClin- 
tock, 2006). Observationally, this is characterized by a 
cool, faint thermal disk and a hard power-law compo- 
nent (photon index F « 1.5-2). The power-law emis- 
sion comes from the radiatively inefRcient inner region 
of the inflow, and is probably due to inverse-Compton- 
scattering in a hot, optically-thin corona or at the base 
of a jet. The disk is sometimes truncated at a large 
distance from the innermost stable circular orbit. Ra- 
dio studies suggest (Fender et al., 2004) that most of 
the accretion power is carried out by a steady jet (de- 
tected as flat-spectrum radio emission), with a possible 
advected component. 

When 0.1 < m < 1, accreting BHs are typically dom- 
inated by an optically-thick, geometrically-thin disk. 



2 



emitting a multicolour blackbody spectrum (Shakura 
and Sunyaev, 1973); this is known as high/soft state, 
or thermal-dominant state (Remillard and McClintock, 
2006). Most of the accretion power is radiated ef- 
ficiently. Advection is not significant, and jets are 
quenched. 

What happens at m > f (very high state, or steep- 
power-law state) is less well known. Few Galactic BHs 
have reached this threshold, and only for short periods 
of time. The X-ray spectrum becomes dominated again 
by a non-thermal component, with the disk contribut- 
ing ~ 10-50% of the total flux. As in the low/hard 
state, the power-law emission may be due to Comp- 
tonization of the disk photons. Radiatively driven out- 
flows are also expected, and radio flaring is observed in 
some sources. It is not clear whether the inner disk is 
covered or replaced by outflows, or by a corona. It is 
also not clear how the accretion power is transferred to 
or directly released in the upscattering medium. 

A strictly related problem is the interpretation of 
ultraluminous X-ray sources (ULXs). Both their X- 
ray spectral and timing behaviour have similarities with 
the very high state. It was also suggested that ULXs 
may have high accretion rates m 3> 1 (Begelman et al., 
2006) or BH masses much higher than those of Galactic 
stellar-mass BHs (Miller et al., 2004). However, a clear 
understanding of the quantitative connection between 
ULXs and stellar-mass BHs is still missing. Here, we 
try to explore this comparison and determine whether 
ULXs are a separate class of BHs, or are stellar-mass 
BHs in a different accretion state. 



2 Thermal disk equations for the high/soft 
state 

We start by assuming that the disk extends all the way 
to the innermost stable circular orbit, so that 



where: Tin is the peak colour temperature; T\^l k is the 
peak effective temperature; k is the hardening factor 
(1.5 ^ K < 2.6); Tin is the "apparent" inner-disk ra- 
dius. The numerical factor ^ was introduced (Kubota 
et al., 1998) to obtain a correctly normalized bolomct- 
ric disk luminosity, taking into account that the fitted 
peak temperature occurs at i? = (49/36)i?in because of 
the no-torque boundary conditions. With this approx- 
imation, and with the choice of k = 1.7 (Shimura and 
Takahara, 1995), the physical inner-disk radius Rin is 
related to the apparent value rin by: 



1.19ri, 



(4) 



The bolometric disk luminosity can be obtained from 
the observed (or extrapolated) flux: 



idisk = 27r(i^/boi(cosz) \ 



(5) 



where i is the viewing angle (i = means face-on). 
can be directly inferred from X-ray spectral fitting^; rir 
is derived from Ldisk and Tin via Equation (3). 
From Eqs. (1), (3), and (4) it follows that: 
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which can also be expressed as Ldisk ~ M^T-^. This is 
the fundamental evolutionary track of an accretion disk 
in the thermal-dominant state. Despite the various ap- 
proximations, Eq. (6) works well (within a factor of 2) 
when applied to the masses of Galactic BHs, assuming 
radiative efficiencies < 0.2. 



Rin = i?isco = 6aGM/c' 



(1) 



where a depends on the spin of the BH (a = 1 for a 
Schwarzschild BH, a = 1/6 for an extreme Kerr BH). 
The bolometric disk luminosity is directly related to the 
accretion rate, from general energy-conservation prin- 
ciples. Ignoring relativistic corrections. 
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where rj is the radiative efficiency. 

In the standard disk-blackbody model (Makishima 
et al., 1986, 2000), the bolometric disk luminosity is 
approximated by 



^disk 



47ra(Tin/K)4(i?in/e)' = 4^'T7l4r?„ 



(3) 



3 Tracks in the luminosity— temperature plane 

3.1 Evolution along the thermal track 

As its accretion rate varies, a source in the high/soft 
state moves along a track of constant M and Rin- This 
track typically extends for about one order of magni- 
tude in luminosity (e.g., Kubota and Makishima, 2004; 
Miller et al., 2004). At the upper end, it is truncated 
where the disk luminosity approaches the classical Ed- 
dington limit, ^Edd ~ 1-3 x 10^^{M/Mq) erg s"i. At 



^the most commonly used implementation of the disk-blackbody 
model is diskbb in XSPEC (Arnaud, 1996). 
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the lower end, it ends as the source reverts to the 
low/hard state. 

Another process that causes the source parameters 
to slide along the thermal track is the removal of a 
fraction of power e > 0, extracted from the disk at 
each radius (independent of radius). For example, this 
energy may be used to power a corona. As a result, 
the peak temperature Tin — > = (1 — e)^/^Ti„, and 
idisk idisk = (1 ^ e)idisk- Neither Rin nor ri„ 
changes; the emitted spectrum is still a disk blackbody. 
This scenario is equivalent to changing the effective ac- 
cretion rate M M' = (1 — e)Af. In summary, changes 
in the accretion rate, or the (uniform) draining of power 
via non-radiative processes, do not affect the BH mass 
estimate from the fitted disk parameters via in Eq. (6). 

3.2 Evolution off the thermal track 

As the accretion rate approaches or exceeds the Edding- 
ton limit, an increasing fraction of the photons emitted 
by the disk are mildly upscattered by hot electrons at 
the disk photosphere. This can be modelled by increas- 
ing the hardening factor k, typically from w 1.7 to w 2.5 
(Davis et al., 2005; 2006). In the (Tin,Ldisk) plane, the 
effect is to move the source horizontally to the right: 
higher colour temperature Ti„ but same effective tem- 
perature Tin/K and same luminosity .Ldisk- The true 
inner-disk radius i?in does not change, while the appar- 
ent disk radius decreases: rin = [k/{k + AK)]'^rin- 
This does not affect the mass estimate, which depends 
on Tin/n. 

Another effect that needs to be taken into account 
is a partial covering of the disk surface: a fraction X 
of photons may be absorbed or more generally removed 
from the disk-blackbody spectrum by clouds or a mod- 
erately optically-thick corona. They are re-emitted in 
other spectral bands or components (e.g.. X-ray power- 
law, or infrared). In a self-consistent analysis, one must 
account for all upscattering, downscattering and ab- 
sorption effects, redistributing the (fixed) total avail- 
able accretion power into the various components. But 
for the purpose of relating the disk emission param- 
eters to the BH mass, we can simply consider those 
photons as being lost from the disk spectrum. The ef- 
fect is to reduce the observed disk flux and luminosity: 
Ldisk —>■ ^disk = (1 ~ -''')-^disk- The spectral shape is 
not altered. The peak colour temperature Tin is not 
changed. Rin stays the same but the apparent disk ra- 
dius decreases: rin r'^^ = (1 ~ Xy/'^rm, because rin 
is indirectly derived from the flux normalization. The 
BH mass inferred from Eq. (6) is a factor (1 — X)^^'^ 
less than the true mass. 

The combined effect of spectral hardening and 
Compton upscattering is to move the location of the 



source to the right-hand-side of its thermal track, in 
the (Tin,Tdisk) plane (Figure 1). Physically, it leads 
to lower estimates for the fitted radius rin (sometimes 
much lower than the innermost stable orbit), and may 
lead to an underestimate of the BH mass if not properly 
accounted for. As expected, accreting BHs occupy this 
region of the parameter space when they are in their 
very high state (or steep-power-law state) (Remillard 
and McClintock, 2006), at accretion rates higher than 
in the high/soft state. 

4 The very high state 

4.1 Very high state in the stellar-mass BH H1743-322 

The Galactic BH H1743-322 provides a textbook ex- 
ample of state transitions (McClintock ct al., 2007). It 
showed a strong outburst between 2003 April and 2003 
October (Markwardt and Swank, 2003; Miller et al., 
2006). Starting from the low/hard state, it remained 
hard (power-law dominated) as it quickly increased in 
luminosity, until it reached the very high state near the 
peak of its outburst (as monitored by RXTE). It then 
oscillated between the very high and the high/soft state, 
eventually settling in the high/soft state near the high- 
luminosity end of the thermal track. Then, it moved 
downwards along the thermal track (decreasing rn) and 
finally returned to the low/hard state (McClintock et 
al., 2007). 

Significantly, the disk parameters were always on the 
right- hand-side of the thermal track in the (Tin, Tdisk) 
plane (Figures 2 and 3), suggesting that we were always 
seeing a full disk, partly modified by Comptonization. 
There is no indication that the inner boundary of the 
disk retreated to larger radii at any stage, either in the 
low/hard or in the very high state (Figure 4), even when 
most of the X-ray flux was carried out by the power- 
law component. The disk appeared subluminous for 
its temperature, compared to the thermal track, be- 
cause some disk photons had been removed from the 
observed thermal component. Spectral fits give higher 
peak colour temperatures, and lower apparent radii, 
when the source was power-law dominated. This be- 
haviour is typical of other (transient) Galactic BHs with 
low-mass donor stars; for example, GRO J1655— 40 dur- 
ing the 1996 outburst (Sobczak et al., 1999). 

4.2 Very high state in ULXs? 

In many respects, ULXs are consistent with being in 
the very high state, based on their X-ray spectral and 
timing properties (e.g.. Done and Kubota, 2006; Dc- 
wangan et al., 2006; Goad et al., 2006; Strohmayer et 



4 



al., 2007). However, one of the reasons it has been 
diffi-Cuh to interpret their nature unequivocaUy is that 
their X-ray spectra can often be fitted equaUy weU, in 
our hmited observing window 0.3-10 keV), with a 
variety of phenomenological or physical models, each 
corresponding to a very different physical scenario. Ex- 
amples of this degeneracy are discussed for example in 
Gongalves and Soria (2006), Stobbart et al. (2006), 
Feng and Kaaret (2007). In addition, sometimes the 
same phenomenological fitting model can have two com- 
pletely different physical interpretations. In the "stan- 
dard" phenomenological model (sometimes known as 
the "cool disk model" ) , which seems to be applicable to 
the majority of sources, ULX spectra have a dominant 
power-law component, with a thermal (disk?) compo- 
nent contributing ^ 10-50% of the X-ray flux (Stobbart 
et al. 2006). Their fitted temperatures are typically a 
few times lower than in stellar-mass BH disks, and their 
apparent inner-disk radii are two orders of magnitude 
higher. Spectral hardening and partial covering of the 
disk cannot explain this difference. If ULXs lie near or 
on the right-hand-side of their thermal tracks, their BH 
masses ought to be two orders of magnitude higher than 
in Galactic systems. This argument has been invoked in 
support of intermediate-mass BHs (Miller et al. , 2004) . 
This physical interpretation of the cool disk model may 
be unsatisfactory, until independent evidence is found 
for the formation of intermediate-mass BHs in the local 
Universe. We need to search for a simpler physical in- 
terpretation of the cool disk model that does not require 
new astrophysical objects. 

5 Ultraluminous branch in XTE J1550-564 

XTE J1550— 564 is one of the best-known microquasars; 
its BH has a dynamical mass « lOM© (Orosz et 
al., 2002). There is some evidence that the source 
was in a low/hard state at the beginning of the 1998 
September - 1999 April outburst, just before the ini- 
tial steep rise (Sobczak et al., 2000). In the early 
phase of the outburst, when the integrated X-ray lu- 
minosity and presumably the accretion rate were at 
their peaks, the X-ray spectrum was dominated by a 
power-law, with a relatively minor disk contribution. 
This, and the presence of characteristic low-frequency 
quasi-periodic-oscillations (LF-QPOs), are typical sig- 
natures of the very high state (Remillard and McClin- 
tock, 2006). However, there was an important differ- 
ence. During the first (brightest) four weeks of the out- 
burst, the source spent some time in the "normal" very 
high state, on the right-hand-side of its thermal track, 
and some time well to the left of that track. This is 
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Fig. 1. — Relation between high/soft state (thermal 
track Ldisk ~ marked by the dashed line) and very 
high state. Source parameters move up along the track 
(A) when m increases; move down (B) when m de- 
creases or when a fraction of power is transferred to an 
optically-thin corona. Partial covering of the disk (C) 
and hardening of the photon spectrum (D) move the 
source parameters to the right-hand-side of the track 
(very high state). From Soria and Kuncic (2007a). 
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Fig. 2. — Spectral evolution of the stellar-mass BH 
H1743— 322 in the temperature-luminosity plane, dur- 
ing the 2003 April-October outburst (data from Mc- 
Clintock et al., 2007). The source parameters are al- 
ways located on the right-hand-side of the thermal 
track; this suggests that the disk extends to i?isco even 
when it is partly modified by a corona. Red datapoints 
mark the thermal-dominant state, as identified by Mc- 
Clintock et al. (2007); green circles are for the very high 
state; blue circles for the low/hard state; black circles 
for intermediate or transition states. We assumed a 
distance of 10 kpc; the viewing angle is unknown. 
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Fig. 3. — Close-up view of the high/soft and very high 
states in the outburst of H1743— 322; colours are de- 
fined as in Figure 2. The dashed line marks the thermal 
track. As the accretion rate increases, the source de- 
viates more and more from the thermal track (see also 
Figure 1). 
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Fig. 4. — Apparent inner-disk radius for H1743— 322, 
plotted as a function of the peak colour temperature 
of the disk. Colours are defined as in Figure 2. In 
the high/soft state, nn « -Risco (Eq- 4). In all other 
states. Tin < ^iscO) suggesting a full disk corrected by 
a variable hardening factor due to Comptonization. 



clearly noticeable both in the original set of spectral 
models (Sobczak et al., 2000, which we reproduce here 
in Figure 5), and in the re-analysis of the data done 
by Kubota and Done (2004) with more complex spec- 
tral models. When XTE J1550-564 was on the left of 
the track, the peak temperature of the disk was much 
cooler than in the normal very high state, and the fit- 
ted inner-disk radius was much larger than the inner- 
most stable orbit. This behaviour has not been seen be- 
fore in other stellar-mass BHs, but may provide a clue 
to understand ULXs. The new state was interpreted 
(Kubota and Done, 2004; Done and Kubota, 2006) as 
a "strong very high state" in which the disk may not 
extend to the innermost stable orbit or may be covered 
by a dominant comptonizing corona. We refer to this 
new state as "ultraluminous branch" (Figures 5, 6). 

At a later stage of the outburst, as the mass accretion 
rate was probably declining, XTE J1550— 564 moved 
back to the right-hand-side of the thermal track, and os- 
cillated a couple of times over the following few months 
between this state and the high/soft state (perhaps in 
response to a moderatley variable but still rather high 
accretion rate) . Eventually, the source declined, moving 
downwards along the thermal track. As it re-entered a 
low/hard state, the disk parameters moved again to the 
left of the thermal track (Figure 5, 6). 

6 Physical interpretation of the ultraluminous 
branch 

Spectral hardening and partial upscattering of the disk 
photons do not explain why the disk should appear 
larger and cooler, at such high X-ray luminosities; one 
might expect the opposite behaviour. Since the BH 
mass does not change, a large apparent radius is gen- 
erally explained as evidence of a truncated or invisible 
inner disk. Indeed, it has been suggested (Done and 
Kubota, 2006; Dewangan et al., 2006; Goad et al., 2006) 
that when the source is on the ultraluminous branch, all 
the photons emitted from the inner disk are upscattered 
by an optically-thick, low-temperature corona. What 
we still see as direct disk emission comes from further 
out, beyond a transition radius Rc » ^iscOj where the 
disk is also cooler. Another possibility is that the in- 
ner disk is disrupted or simply occulted by an optically- 
thick, radiatively-driven outflow (Poutanen et al., 2006; 
Begelman et al., 2007). Alternatively, the inner disk 
may still be present and directly visible, but most of 
its accretion power is extracted via non-radiative pro- 
cesses, such as mechanical outflows or Foynting flux 
(Kuncic and Bicknell, 2004). In this case, the peak 
temperature occurs much further away from the inner 
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Fig. 5. — Luminosity-temperature plot for the stellar- 
mass BH XTE J1550-564 during its 1998 September - 
1999 April outburst. Accretion states are colour-coded 
as in Figure 2. What was traditionally classified as the 
very high state (green datapoints) is split here into two 
states, approximately defined as being on the left (ul- 
traluminous branch) and on the right (very high state) 
of the thermal track. Arrows indicate the approxime 
evolution track during the outburst. The fit parame- 
ters are from Sobczak et al. (2000), and we assumed a 
distance of 5 kpc and a viewing angle i ~ 70°. 




kTi„ (keV) 

Fig. 6.— Inner-disk radius of XTE J1550-564 plot- 
ted as a function of peak colour temperature. Accre- 
tion states are colour-coded as in Figure 2. A constant 
conversion factor Rm = 1.19rin was used for all dat- 
apoints. The radius becomes :s> i?isco ~ 60-70 km 
when the source is on the ultraluminous branch; the 
disk temperature decreases along the same branch. An 
apparently retreating disk is also found in the low/hard 
state. 



disk edge. Some of the power released through non- 
radiative channels would then be converted to power- 
law photons. We leave a more detailed discussion of 
the physical interpretation of the transition radius to 
further work (Soria and Kuncic, 2007b). 

An important question to ask is whether the ultra- 
luminous branch corresponds to a higher or lower ac- 
cretion state than the classical very high state. The de- 
tected 2-20 keV flux {RXTE da.ta., Sobczak et al., 2000) 
was higher in the very high state. However, that is 
not necessarily a good indicator of the total luminos- 
ity, especially when the disk emission is very soft, nor 
of the accretion rate. For example, the total luminos- 
ity at the top of the thermal track is w ^Edd, with a 
radiative efficiency > 0.1. The non-thermal processes 
responsible for the dominant power-law emission in the 
very high and ultraluminous states are expected to have 
lower radiative efficiency. So, the total luminosity may 
still be w ^Edd, just differently redistributed between 
thermal and non-thermal components, even if m > a 
few. In the case of XTE J1550— 564, the total disk plus 
power-law flux is at least as high, and probably slightly 
higher, along the ultraluminous branch (Soria and Kun- 
cic, 2007b). More importantly, the disk luminosity is 
slightly increasing along the ultraluminous branch, for 
an increasing inner-disk radius. This can happen (recall 
Eq. 2) only if the accretion rates also increases along 
the ultraluminous branch (even if the total luminosity 
did not). The fact that XTE J1550-564 reached both 
the very high state and the ultraluminous branch in 
the first, more energetic phase of its outburst, and only 
(briefiy) the very high state in the second, weaker phase 
(Sobczak et al., 2000) also suggests that the ultralumi- 
nous branch corresponds to a higher activity state. As 
the source started to decline, it went from very high 
state to high/soft state, and finally to low/hard state, 
consistent with this interpretation. 



7 Disk parameters on the ultraluminous 
branch 

A plausible phenomenological way to model the disk 
parameters on the ultraluminous branch is to assume a 
transition radius i?c S> i?isco- For R > Rc, the inflow 
can be approximated by a standard disk, such that 



L. 



disk 



47rcrTi'i?2 



(7) 



where Tc = T{Rc) is the maximum (observable) disk 
temperature. For i?isco < i? < -Rc, we assume that 
all the accretion power comes out in the power-law-like 
component, and the standard disk either is not directly 
visible, or is disrupted, or emits a negligible direct flux. 



Bridging the gap between stellar-mass black holes and ultraluminous X-ray sources 



7 



So, regardless of the details, Rc will appear as the inner 
radius when we fit the disk spectrum. In this scenario, 
Eq. (1) no longer holds or is no longer relevant to the 
observed spectrum. It is replaced by 



Rc = FRisco 



(8) 



with F = F{m) > 1. The visible disk radiates only 
the accretion power released from the outer radius to 
Rc- The effective radiative efficiency of that part of the 
inflow - M/Rc ^ l/(12aF) < 0.1. Therefore, we have 
(cf. Eqs. 2 and 6): 
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where the first expression for Af is scaled to the typ- 
ical values observed in stellar-mass BHs, and the sec- 
ond one is scaled to typical ULX parameters. An esti- 
mate of the BH mass based only on the disk luminosity 
and temperature, without taking into account the trun- 
cation factor F, will necessarily over-estimate the BH 
mass. We argue that this is the main reason why the 
intermediate-mass BH hypothesis is probably incorrect. 

Hence, to estimate the BH mass and accretion rate of 
a source on the ultraluminous branch (Eqs. 9 and 10), 
we need to measure F — Rcl Risco directly, or model 
its dependence on rh. For XTE J1550— 564, we see that 
the apparent radius increases by a factor F « 6-8 along 
that branch (Figure 6). For ULXs, we do not have a 
complete coverage of state transitions, and we do not 
know the BH mass independently. The simplest way to 
estimate F, at least as an order of magnitude, is to as- 
sume that all of the accretion power released at R> Rc 
is radiated by the disk, and some or most of the power 
released at -Rjsco < < i?c is eventually radiated as 
power-law photons, without contributing to the ther- 
mal disk component. From Eqs. (2) and (9), we get 
F > itot/^disk- The inequality sign takes into account 
the fact that the radiative efficiency of the non-thermal 
processes in the inner region is less than the efficiency 
of blackbody emission in a standard disk. For the most 
luminous nearby ULXs, extrapolating from the relative 
fluxes of the power-law and thermal components in the 
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Fig. 7.— As in Figure 5, but with two ULXs (NGC 1313 
X-1 and X-2) plotted alongside the stellar-mass BH 
XTE J1550— 564. Both ULXs show an anticorrelation 
between disk luminosity and temperature. This sug- 
gests that they are also on their respective ultralumi- 
nous branches. Unlike XTE J1550— 564, they never oc- 
cupy any other accretion state. The colour coding of 
the XTE J1550— 564 datapoints is the same as in Fig- 
ure 2. Datapoints (and, where available, error bars) 
for NGC 1313 X-1 are plotted in magenta; squares are 
for the XMM-Newton data (Feng and Kaaret, 2006, 
2007), and diamonds for the Suzaku data (Mizuno et 
al., 2007). Datapoints for NGC 1313 X-2 are plotted in 
orange (triangles for the XMM-Newton data and a star 
for the Suzaku data). 

0.3-10 keV band, the luminosity ratio Ltot/^disk ~ 3- 
20; in a few other, more distant ULXs, there is only a 
lower limit > 10 for this ratio (Stobbart et al., 2006; 
Winter et al., 2006). Hence, we infer that F > 10 for 
typical ULXs. 



8 Observational predictions 

From the fitted disk luminosities and temperatures 
(e.g., Miller et al., 2004; Feng and Kaaret, 2005; Stob- 
bart et al., 2006), typical values of i?c > 5000 km are 
found for many ULXs, i.e., ~ 100 times larger than typ- 
ical inner-disk radii in stellar-mass BH. If, as we spec- 
ulate, ULXs are accreting BHs on the ultraluminous 
branch, with F > 10, it means that their BH masses 
are only required to be < 10 times larger than typical 
BH masses in Galactic systems. This simple argument 
suggests that ULX masses are < IOOMq. At the same 
time, masses ~ 5O-1OOM0 are still sufficiently high to 
guarantee that the apparent ULX luminosity is com- 
parable or only a factor of a few higher than the Ed- 
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Fig. 8.— As in Figure 6, but with two ULXs (NGC 1313 
X-1 and X-2) plotted alongside the stellar-mass BH 
XTE J1550— 564. A constant conversion factor Ri^ = 
1.19rin was used for all datapoints. A viewing angle 
i = 60° was assumed for NGC 1313 X-1 and X-2, for 
consistency with the flux/luminosity conversion of Feng 
and Kaaret (2006, 2007). The colour coding is the same 
as in Figures 7. The anticorrelation of radius and tem- 
perature seen in the two ULXs is very similar to (in 
fact, almost the direct extension of) the ultraluminous 
track of XTE J1550-564. 

dington limit, and therefore they do not require strong 
collimation of their emission. 

The location or evolutionary track of a source in the 
(Tci^disk) plane provides another observational con- 
straint. To predict the displacement of a source from its 
thermal track when F > 1, we need to understand how 
F, Tc and Ldisk vary as a function of to. For a fixed BH 
mass, the radiative flux equation of a standard disk tells 
us that Tc ~ i?^r^^^TO^^^- It was suggested (Poutanen 
et al., 2007; Begelman et al., 2007) that Rc ^ m, based 
on plausible physical processes that may form such a 
transition radius. If so, from Eq. (7) we expect that, as 
the accretion rate m increases, a source will move along 
a track parameterized by Ldisk ~ constant, Tc ^ m~^/^. 
This is not exactly what is found in the RXTE spec- 
tra of XTE J1550-564, because, there, Ldisk - T-^- 
However, the luminosity track is very sensitive to the 
radial temperature distribution on the disk, at R > Rc- 
Disk models with a distribution flatter than are 
sometimes used. For example, it was found (Kubota et 
al., 2005) that T i?""-^ may provide a more accu- 
rate fit to the X-ray spectral data of Galactic BHs. In 
that case, we expect Ldisk j'-o.44^ with Tc ~ m~°-^^. 
For a disk dominated by X-ray irradiation (as may be 
the case it R ^ ^isco), T ^ R~°-^. In that case, 
-C'disk T^^, with Tc ~ TO~°-^^. Another way to re- 



produce the observed anticorrelation of disk luminosity 
and peak temperature is to parameterize Rc ^ with 
/? < 1. For example, Ldisk ~ T^^-^ for Rc ^ w?!^ and 
a standard temperature profile. 

Regardless of the precise functional forms of F[jn) 
and Tc(i?c), the significant result is that we expect a 
source to move to the left-hand- side of its thermal track, 
as the accretion rate increases to to ^ 1. Along that 
track, the luminosity of the disk component stays con- 
stant or increases, even if the observed peak tempera- 
ture decreases. Thus, it appears that in order of increas- 
ing accretion rates, the spectral states of an accreting 
BH are: high/soft (thermal track); very high (to the 
right of that track); ultraluminous (to its left). Along 
this sequence, we also expect the source to become more 
and more dominated by non-thermal emission compo- 
nents. 



9 Testing the ultraluminous scenario vifith 
fitted ULX parameters 

Let us consider an accreting BH with M « 5OM0. 
If TO « 1 and F — 1, we expect its spectrum to be 
dominated by the disk component, with a peak tem- 
perature ~ 0.8 keV, a characteristic size of the X-ray 
emitting region ^ GGM/c^ ~ 500 km, and a disk lumi- 
nosity « 5 X 10'^^ erg s^^ « ^Edd- This luminosity is 
consistent with those found in average ULXs, but the 
other fit parameters are not. Let us now suppose that 
T « TO « 20, for the same source on the ultraluminous 
branch. In that case, we expect that the source will still 
have idisk ~ ^Edd « 5 X 10'^^ erg s~^ or shghtly higher, 
but a peak temperature « 0.8 x 20^°-^ keV « 0.18 
keV, a characteristic radius ^ 20 x 6GM/c^ ~ 10'' km, 
and an X-ray power-law component with a luminosity 
> 5 X 10^^ erg s~^. These are more typical ULX pa- 
rameters. If we had observed such a source, and had 
directly inserted the fitted values of Tdisk and Tc into 
Eq. (6) instead of Eq. (10), we would have incorrectly 
interpreted the source as an intermediate-mass BH with 
Al « 1OOOM0. The same mistake would have occurred 
if we had only observed XTE J1550— 564 at the end of 
its ultraluminous branch. In that case, we might have 
estimated a mass « 60-80Mq from Eq. (6), instead of 
1OM0. 

We want to test whether at least some ULXs show 
evidence of spectral evolution along the ultralumi- 
nous branch, with a predicted (and observed in XTE 
J1550— 564) anticorrelation between disk luminosity 
and temperature. For this, we consider the two "clas- 
sical" sources, NGC 1313 X-1 and X-2; XMM-Newton 
and Suzaku observations of those sources were recently 
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studied by Feng and Kaaret (2006, 2007) and Mizuno 
et al. (2007). Their works show that both sources ad- 
mit at least two alternative fitting models: a cool-disk 
and a hot-disk scenario. The hot-disk model was con- 
sidered more physical, partly because it would place 
those ULXs close to the extrapolation of a thermal track 
(-^disk ~ while the cool-disk model place them to 
its left and would imply an anticorrelation between disk 
temperature and luminositjU. However, based on the 
comparison with XTE J1550— 564 (Figure 7), this is 
precisely what we expect to find. Therefore, we argue 
that the spectral evolution of those two ULXs strongly 
supports the cool-disk model. The physical interpre- 
tation of the cool disk is not based on intermediate- 
mass BHs, but on a receding inner radius i?c, when the 
sources are moving along their ultraluminous branches 
at accretion rates m > 10. The increase in the appar- 
ent inner-disk radius Rc is even more obvious when it 
is plotted against Tc (Figure 8). We suggest that this 
is further evidence of a fundamental similarity between 
the ultraluminous branch seen in XTE J1550— 564 and 
the ULX behaviour. Note that, unlike XTE J1550-564, 
neither NGC 1313 X-1 nor X-2 have shown evidence of 
ever reaching their thermal tracks (Figures 7, 8), which 
remain undetermined. This makes it more difficult to 
disentangle the effects of higher masses and higher ac- 
cretion rates, and hence to estimate the values of F and 
of their BH masses. We will discuss these issues quan- 
titatively and in greater details in a forthcoming paper 
(Soria and Kuncic, 2007b). 

10 Evidence from X-ray timing: low-frequency 
QPOs 

Some nearby, luminous ULXs show LF-QPOs in their 
power-density-spectra (Strohmayer et al., 2007; Dewan- 
gan et al., 2006; Strohmayer and Mushotzky, 2003). 



caveat about the cool-disk fitting model for NGC 1313 X-1 
and X-2 is that the disk temperature also appears to be anticor- 
related with the column density (see in particular Fig. 1 of Feng 
and Kaaret 2007). Thus, it has been suggested that the apparent 
increase in the disk luminosity at lower disk temperatures could 
be a spurious effect, if the anticorrelation between column den- 
sity and disk temperature has no physical basis. However, we 
argue that an increase in the intrinsic absorption as the accre- 
tion rates increases to more than one order of magnitude above 
Eddington, along the ultraluminous branch, is in fact very plau- 
sible. Other transient X-ray binaries show higher absorption near 
their outburst peaks. So, we do not think this is a serious reason 
to consider the disk luminosity-temperature anticorrelation un- 
physical. It is, however, possible that the luminosity increase at 
lower temperatures is partly due to uncertainties in the column 
density, and therefore that the slope of the ultraluminous branch 
is flatter than it appears in Figure 7. 
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Fig. 9. — Frequencies of the LF-QPOs detected (only) 
in the very high and ultraluminous states of XTE 
J1550-564 during its 1998 September - 1999 April out- 
burst, as a function of fitted inner-disk radius. A con- 
stant conversion factor i?in ~ 1.19rin was used for all 
datapoints. The anticorrelation of frequency and radius 
(and hence, of frequency and accretion rate) provides 
an independent constraint to ULX masses and their ac- 
cretion rates. 

Such oscillations are more strongly associated with the 
power-law photons (Strohmayer et al., 2007). A simi- 
lar behaviour is observed in stellar-mass BHs, includ- 
ing XTE J1550— 564, in their very high and ultralumi- 
nous states (Remillard and McClintock, 2006). This 
has been interpreted (e.g.. Done and Kubota, 2006) as 
further evidence of a connection between those lumi- 
nous, power-law dominated states in ULXs and stellar- 
mass BHs, also consistent with their spectral properties. 
Moreover, it was suggested that LF-QPO frequencies 
scale as 1 /M across the whole range of BH masses, from 
stellar-mass objects to AGN (McHardy et al., 2006). 

LF-QPO frequencies in ULXs are a factor « 20-100 
lower than in stellar-mass BHs; for example, one of the 
most notable cases is the 20-mHz QPO in NGC 5488 
X-1 (Strohmayer et al., 2007). At face value, this 
appears to strengthen the intermediate-mass BH hy- 
pothesis (Strohmayer et al., 2007). However, just as 
we discussed for the large inner-disk radii, there may 
be a more physical interpretation, which becomes ap- 
parent from a comparison with the LF-QPOs in XTE 
J1550— 564. From the timing analysis of Sobczak et al 
(2000), it is clear (Figure 9) that the LF-QPO frequency 
in this stellar-mass source decreased as the apparent ra- 
dius increased, during the 1998-1999 outburst. That is, 
for a fixed BH mass, the frequency decreased as the ac- 
cretion rate increased along the proposed ultraluminous 
branch. The relation is not linear; however, it can be 
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roughly approximated as i^qpo ~ rh""'^ over most of 
the parameter range. More generally, taking into ac- 
count the expected (1/M) scaling between sources of 
different masses, we can write 

i^QPO ~ rh-^ (1/M) (11) 

with 0.5 < (5 < 1. 

This suggests that the lower frequencies found in 
ULXs may also be the result of higher accretion rates, 
not only of higher masses, in agreement with our pro- 
posed interpretation of the spectral data. The different 
dependence of spectral and timing parameters on m and 
M suggests that it is in principle possible to disentan- 
gle the two effects. Preliminary back-of-the-envelope 
calculations suggest that accretion rates m « 20 and 
masses M ~ SO-IOOM© are consistent with both the 
spectral and timing data of typical ULXs. 

In addition to providing an independent constraint 
on the effects of accretion rate and BH mass, the ob- 
served anticorrelation between a disk parameter (inner 
radius) and the frequency of oscillation of the power-law 
emission component is also a fundamental clue to un- 
derstand the origin of such oscillations. One can spec- 
ulate that LF-QPOs originate at the interface of outer 
disk and inner flow, at the transition radius Re- We will 
present a more detailed analysis of the timing results, 
in the framework of our ultraluminous branch model, 
in Soria and Kuncic (2007b). 

11 Conclusions 

It is often noted that ULXs and stellar-mass BHs share 
X-ray spectral and timing properties typical of the 
luminous, power-law dominated accretion state (very 
high state). If we plot the spectral evolution of typ- 
ical accreting BHs in the (Tin,idisk) plane, the very 
high state is located on the right-hand-side (higher 
disk temperatures) of the thermal track characteristic 
of the disk- dominated high/soft state. This would im- 
ply that ULXs and stellar-mass BHs are entirely sepa- 
rate species, with a two-order-of-magnitude gap in BH 
masses between them (intermediate-mass BH scenario). 
However, RXTE spectral fits to the stellar-mass BH 
XTE J1550-564 in its 1998-1999 outburst (Sobczak et 
al. 2000; Kubota and Done 2004) clearly show that the 
very high state can in fact be divided into two sub- 
states, located to the right (higher temperature) and to 
the left (lower temperature) of the thermal track, re- 
spectively. In the latter state (ultraluminous branch), 
the fitted inner-disk radius is much larger than the in- 
nermost stable circular orbit, and the peak colour tem- 
perature much cooler than in the high/soft state. This 



is consistent with a standard outer disk truncated or 
obscured beyond a transition radius Rc » Risco- The 
inner region of the inflow contributes mostly to the 
power-law component, perhaps through upscattering in 
a moderately optically-thick corona, or in a magnetized 
wind (Kuncic and Bicknell, 2004). The observed an- 
ticorrelation between bolometric disk luminosity and 
peak colour temperature, and the observed sequence of 
spectral state transitions imply that the ultraluminous 
branch occurs at even higher accretion rates than the 
classical very high state, hence probably at M > a few 

MEdd- 

The ultraluminous branch provides a bridge between 
stellar-mass BHs and ULXs. It supports the cool-disk 
spectral model for ULXs, but not its interpretation in 
terms of intermediate-mass BHs. At the other end 
of the bridge, we showed that two prototypical ULXs 
(NGC 1313 X-1 and X-2) also move along a track con- 
sistent with the ultraluminous branch, with anticorre- 
lations between disk luminosity and temperature, and 
between inner-disk radius and temperature. If this in- 
terpretation is correct, all the main spectral and timing 
properties of accreting BHs (including their characteris- 
tic LF-QPO frequencies) depend simultaneously on two 
factors: a mass scaling and an accretion rate scaling, at 
least when M > a few Msdd- We suggest that the 
observed ULX parameters may be best explained with 
BH masses ~ 50-100Mq and M « 20MEdd- This also 
implies that the radiative efficiency of the non-thermal 
medium is < 0.02. 

In this scenario, ULXs persistently occupy the high- 
accretion-rate ultraluminous branch, whereas stellar- 
mass BHs would only rarely reach such rates, and 
only near the peak of their outbursts (for example, 
we showed in Section 4 that the stellar-mass BH 
H1743— 322 never does). We are currently investigating 
whether brief transitions to an ultraluminous branch 
may have occurred in other historical outbursts of 
Galactic BHs; and whether other variable ULXs move 
along similar tracks. A significant difference between 
the two classes of objects may be caused by the different 
types of Roche-lobe mass transfer. ULXs may accrete 
a few Mq over < 10^ yr from B-type donor stars, while 
transient Galactic BHs accrete from older solar-mass 
stars, which do not persistently fill their Roche lobes. 
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